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Introduction 

A  research  program  consisting  of  basic  experimental  studies  of  the  fluid  mechanics  of 
an  accelerating  airfoil  undergoing  a  stepwise  variation  in  the  angle  of  attack  was  initiated 
in  July  1986  under  grant  no.  AFOSR-86-0243.  This  technical  report  is  being  written  to 
describe  the  work  carried  out  during  the  period,  July  1986  -  December  1987  of  the  grant. 

The  present  study  has  the  following  main  objectives  to  the  general  problem  of  the 
unsteady  vortical  flows  generated  by  high  angle  of  attack  airfoils. 

1.  Investigate  experimentally  the  unsteady  flow  structure  generated  by  a  NACA  airfoil 
at  different  angles  of  attack  started  impulsively  form  rest. 

2.  Investigate  experimentally  the  transient  flow  characteristics  of  an  NACA  airfoil  un¬ 
dergoing  a  stepwise  varied  angle  of  attack. 

3.  Investigate  the  unsteady  flow  structure  of  an  accelerating  NACA  0012  airfoil  at  dif¬ 
ferent  fixed  angles  of  attack. 

4.  Investigate  the  transient  flow  characteristics  of  an  accelerating  airfoil  undergoing  a 
stepwise  incidence  variation. 

As  originally  envisaged  in  all  these  studies,  attention  was  directed  to  the  basic  under¬ 
standing  of  the  unsteady  flow  phenomena.  During  these  investigations,  a  paralleled  effort 
was  devoted  in  consultation  with  AFOSR,  to  study  some  of  the  experimentally  observed 
features  using  numerical  solutions.  A  unique  experimental  technique,  developed  in  our  lab¬ 
oratory,  is  successfully  implemented  to  study,  in  detail,  the  unsteady  large  scale  vortical 
motions  that  occur  in  these  flows. 

In  the  following,  the  status  of  the  research  effort  is  given  along  with  some  of  the  most 
important  conclusions  arrived  in  this  study  to  date. 
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Status  of  the  Research  Effort 

In  the  following,  major  conclusions  on  the  various  facets  of  the  research  program  are 
given. 

A  systematic  investigation  has  been  carried  out  on  the  structure  of  an  unsteady  flow 
field  generated  by  an  impulsively  started  NACA  0012  airfoil,,  of  finite  aspect  ratio,  at 
different  angles  of  attack  (0  <  a  <  45°),  and  at  a  fixed  Reynolds  number  of  1400. 

A  novel  experimental  technique  known  as  “Particle  Image  Displacement  Velocimetry” 
was  used  to  measure  the  instantaneous  two  dimensional  velocity  field.  The  velocity  field 
was  measured  with  sufficient  accuracy  and  spatial  resolution  that  the  vorticity  field  and 
the  pressure  field  can  be  computed  very  accurately,  a  unique  capability  of  this  technique. 
The  detailed  description  of  this  technique  is  given  in  Appendix  I. 

A  parallel  computational  study  was  conducted  to  augment  the  above  mentioned  ex¬ 
perimental  studies.  This  effort  was  supported  by  the  Mechanical  Engineering  Department. 
In  this  study,  the  performance  of  the  discrete  vortex,  random  walk  approximation  to  the 
Navier-Stokes  equations  was  examined  for  the  case  of  large  scale  separation  flow  about  an 
impulsively  started  NACA  0012  airfoil  at  different  angles  of  attack.  The  primary  objective 
was  to  include  sufficient  vortices  to  resolve  the  shortest  meaningful  length  scales  of  the 
flow  at  a  given  computational  time  step.  The  new  fast  velocity  summation  algorithm* 
enables  the  flow  to  be  computed  with  much  more  resolution  than  previously  possible  in 
vortex  methods. 

The  main  features  of  the  unsteady  large  scale  separated  flow  about  an  impulsively 
started  airfoil  are  as  follows: 

The  multiple  exposed  photographs  of  the  flow  field  about  the  airfoil  at  10°  or  less 
incidence  showed  that  the  flow  is  well  behaved  and  attached  to  the  airfoil  over  the  entire 
period  of  observation.  However,  at  large  angles  of  attack  a  >  20°,  the  flow  separates  on 
the  upper  surface  of  the  airfoil  and  generates  large  scale  vortices.  The  following  scenario 
develops  in  time  on  the  upper  surface  of  the  airfoil.  At  the  start  of  the  airfoil,  a  vortex,  at 
the  trailing  edge,  commonly  designated  as  “starting  vortex”  is  generated  and  carried  away 
from  the  airfoil.  Concomitant  with  this  is  the  generation  of  a  separation  bubble  at  the 
leading  edge  of  the  airfoil.  At  a  later  time,  the  separation  bubble  grows  into  an  isolated 

*  van  Dommelen,  Leon,  and  Rundensteiner.  Elke  A.,  “Fast,  Adaptive  Summation  of 
Point  Forces  in  the  Two-Dimensional  Poisson  Equation.”  Submitted  to  the  Journal  of 
Computational  Physics. 
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primary  vortex  with  secondary  vortices  following  behind  it.  This  multiple  vortex  structure 
continues  to  grow  together  and  move  along  the  upper  surface  until  it  reaches  the  trailing 
edge.  At  this  point,  the  primary  vortex  induces  a  vortex  at  the  trailing  edge.  Finally,  the 
primary  vortex  and  the  induced  trailing  edge  vortex  interact  and  generate  a  complex  flow 
field.  However,  for  finite  aspect  ratio  airfoils  or  wings,  a  different  type  of  flow  field  seems 
to  emerge  at  later  stages  of  development.  The  various  events  described  above  occur  at 
different  times,  depending  upon  the  angle  of  attack  and  Reynolds  number. 

Typical  PIDV  measurements  of  the  instantaneous  velocity  field,  at  different  times,  for 
an  airfoil  a  =  30°  are  shown  in  figure  1.  The  aspect  ratio  of  the  wing  was  about  3.  The 
airfoil  travels  from  right  to  left.  The  data  is  presented  in  a  body  fixed  reference  frame. 
The  length  of  the  velocity  vector  corresponds  to  its  magnitude.  The  dimensionless  time  t* 
is  defined  as  where  U  is  the  free  stream  velocity  and  C  is  the  airfoil  chord.  The  starting 
vortex  (at  the  far  right  of  the  picture)  and  the  initial  separation  bubble  at  the  leading  edge 
can  be  seen  clearly  in  the  picture  corresponding  to  t*  =  0.68.  The  primary  vortex  with 
secondary  vortices  following  behind  it  can  be  seen  in  the  figure  at  <*  =  2.02.  The  trailing 
edge  vortex  can  be  seen  in  the  figure  at  t*  =  3.02.  At  t*  >  3,  the  primary  vortex  abruptly 
moves  away  from  the  upper  surface  leaving  behind  a  “vortex  sheet-”  like  structure.  Such  a 
behavior  is  attributed  to  the  inference  of  tip  vortices  which  are  generated  due  to  the  finite 
aspect  ratio  of  the  airfoil.  At  later  times,  for  example  at  t*  —  4.85,  the  tip  vortices  interact 
with  the  separated  flow  on  the  upper  surface  and  generate  a  complicated  three-dimensional 
flow  field.  The  nature  of  this  interaction  and  the  parameters  that  govern  such  a  flow  field 
is  not  known  yet. 

Typical  two-dimensional  computational  results  from  random-walk  vortex  simulations 
of  the  full  Navier-Stokes  equations  are  shown  in  figure  2.  The  angle  of  attack  and  the 
Reynolds  number  are  the  same  as  those  in  the  experiment;  the  results  of  which  are  shown 
in  figure  1.  The  streamline  pattern,  along  with  vorticity,  which  is  represented  in  bit-mapped 
graphics  as  half  tones  are  shown  in  the  figure.  Except  for  the  effect  of  the  finite  aspect 
ratio  of  the  airfoil,,  the  stream  line  pattern  looks  very  much  similar  to  those  observed  in 
figure  1.  To  further  evaluate  these  results,  the  locus  of  the  primary  vortex  as  its  develops 
in  time  is  shown  in  figure  3.  The  computational  results  agree  well  with  the  experiment 
for  t*  <  2.  Beyond  t *  =  2,  it  is  expected  that  the  experimental  flow  field  was  influenced 
by  the  tip  vortices  making  it  to  be  highly  three-dimensional.  The  coefficients  of  lift  and 
drag  as  obtained  from  the  computations  are  shown  in  figure  4.  As  expected,  the  coefficient 
of  lift  increases  with  t *  up  to  a  point  where  the  primary  vortex  is  attached  to  the  upper 
surface.  For  later  times,  where  the  primary  vortex  leaves  the  upper  surface,  the  coefficient 


of  lift  drops  significantly. 

Results  similar  to  these  described  above  were  obtained  at  a  =  10°,  20°  and  45°.  The 
details  of  this  work  will  be  given  in  a  paper  to  be  submitted  to  the  AIAA  Journal. 
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Abst  ract 

A  novel  experimental  technique, 
commonly  referred  to  as  Laser  Speckle 
Velocimetry  (LSV)  or  Particle  Image 
Displacement  Velocimetry  (PIDV),  is  devel¬ 
oped  for  the  measurement  of  instantaneous 
velocity  fields  in  unsteady  and  steady 
flows.  The  main  advantage  of  this  tech¬ 
nique  is  that  the  velocity  field  is 
measured  with  sufficient  accuracy  so  that 
the  distribution  of  vorticity  can  be  cal¬ 
culated  with  accuracy. 

The  PIDV  technique,  which  is  ideally 
suited  for  the  study  of  unsteady  separated 
flows,  has  been  utilized  to  measure  the 
development  of  the  separated  flow  field 
generated  by  an  high  angle  of  attack 
fa=30°)  NACA  0012  airfoil  started  impul¬ 
sively  from  rest. 

1  .  Introducti on 

For  the  solution  of  many  problems  that 
occur  in  high  angle  of  attack  aero¬ 
dynamics,  it  is  a  necessary  to  have  a 
thorough  understanding  of  the  behavior  of 
unsteady  separated  flows.  Although  much 
progress  has  been  made  in  predicting  the 
steady  flow  phenomenon  with  the  use  of 
numerical  methods,  it  is  still  difficult 
to  predict  unsteady  flows  which  contain 
flow  separation.  The  difficulty  mainly 
arises  from  the  fact  that  these  flows  are 
extremely  complex  and  are  not  amenable  to 
standard  experimental  and  numerical 
techniques.  In  view  of  this,  a  novel 
experimental  technique  is  being  developed 
for  the  measurement  of  instantaneous 
velocity  fields  in  unsteady  and  steady 
fluid  flows.  This  paper  provides  a 
description  of  this  technique  along  with 
its  successful  application  to  the  study  of 
an  unsteady  separated  flow  generated  by  ar 
high  angle  of  attack  airfoil. 

In  an  unsteady  flow,  a  single  photo¬ 
graph  of  the  flow  pattern  at  a  given 
instant  does  not  generally  provide  any 
meaningful  information.  In  order  to 
understand  the  unsteady  flow  phenomena,  it 
is  necessary  to  obtain  both  spatial  and 
temporal  information  of  the  entire  flow 
field.  With  this  in  mind,  optical 
techniques  have  been  widely  used  to 
observe  and  measure  properties  of  flow 
fields  such  as  velocities  and  densities. 
Many  of  these  techniques  are  qualitative 
in  nature,  but  of  great  value  in  guiding 
intuition  and  in  suggesting  ways  to 
investigate  the  problem  by  quantitative 
means.  An  admirable  review  of  the  modern 
optical  techniques  in  fluid  mechanics  is 

’  Associate  pTofessot  and  member  of  AIAA 

Assistant  Professor  and  member  cf  AIAA 


given  by  Lauterborn  and  Vogel1 ,  and  the 
reader  is  referred  to  this  article  for 
details.  A  quan t i ta t i ve  flow  visualiza¬ 
tion  technique  would  be  very  helpful  in 
file  study  of  these  flow  fields.  Attempts 
to  accomplish  this  task  by  tracing  ttie 
streaklines  of  injected  particles41 
usually  can  not  provide  the  spatial 
resolution  that  is  required  and  a  large 
amount  of  labor  is  necessary  to  reduce  the 
data  . 

In  unsteady  separated  flows,  as  will 
be  shown  later,  it  is  often  desirable  to 
obtain  the  vorticity  field,  in  addition  to 
the  velocity  field.  However,  measurement 
of  the  vorticity  exceed  the  present  exper¬ 
imental  capability.  This  difficulty 
arises  from  the  fact  that  vorticity  is  a 
quantity  defined  in  terms  of  local  veloc¬ 
ity  gradients.  In  contrast,  the  currently 
available  measurement  techniques,  such  as 
hot-wire  anemometry  or  laser  velocimetry, 
are  sensitive  only  to  the  local  velocity. 
Hence,  measurements  must  be  made  over 
several  points  and  the  resulting  velocity 
components  are  then  analyzed  by  finite 
difference  schemes.  However,  the  errors 
produced  lay  the  necessary  differentiations 
limit  the  accuracy  and  spectral  range.  In 
addition,  the  spatial  resolution  of  this 
method  is  often  not  sufficient  to  measure 
small-scale  fluid  motions  or  rapidly 
changing  velocity  gradients.  As  a  conse¬ 
quence,  the  measured  vorticity  field  is  a 
type  of  spatially  averaged  estimate  of  the 
actual  vorticity  field.  Finally,  this 
method  provides  information  at  only  a 
single  point.  If  information  on  the 
entire  flow  field  is  required,  measure¬ 
ments  must  be  carried  out  sequentially  one 
point  at  a  time.  This  sequential  method, 
although  laborious,  is  straightforward  in 
applications  involving  steady  flows. 
However,  the  method  becomes  very  diffi¬ 
cult,  if  not  impossible  to  implement,  when 
studying  unsteady  flows.  Direct  measure¬ 
ment  of  vorticity  has  been  tried,  for 
instance,  by  injection  of  spherical 
particles  which  rotate  in  the  flow  with  an 
angular  velocity  proportional  to  the  local 
vorticity4  .  Sucli  methods  suffer  the  same 
drawback  of  insufficient  spatial  resolu¬ 
tion  just  mentioned  and  also  can  be  quite 
compl ex . 

Recently,  a  novel  velocity  measurement 
technique,  commonly  referred  to  as  Laser 
Speckle  Velocimetry  (LSV)  or  Particle 
Image  Displacement  Velocimetry  (PIDV)  has 
become  available.  This  technique  provides 
the  simultaneous  visualization  of  the  two- 
dimensional  streamline  pattern  in  unsteady 
flows  as  well  as  tire  quantification  of  the 
velocity  field  over  an  entire  plane.  The 
advantage  of  this  technique  is  that  the 
velocity  field  can  be  measured  over  an 


entire  plane  of  the  flow  field  simul¬ 
taneously,  with  accuracy  and  spatial 
resolution.  From  this  the  instantaneous 
vorticity  field  can  be  easily  obtained. 
This  constitutes  a  great  asset  for  the 
study  of  a  variety  of  flows  that  evolve 
stochastically  in  both  space  and  time, 
such  as  the  unsteady  vortical  flows  that 
occur  in  rotor  craft  and  high -sing le-of- 
attack  aerodynamics.  For  the  background 
of  this  present  technique,  the  reader  is 
referred  to  reference  5. 

The  principle  of  the  technique  is 
given  in  the  next  section  followed  by  the 
description  of  the  apparatus,  instrumen¬ 
tation  and  procedures.  Section  4  provide 
the  results  and  their  description  for  the 
flow  over  an  high  angle  of  attack  NACA 
0012  airfoil. 


2.  Principle  of  the  Technique 

The  application  of  LSV  or  PIDV  to 
fluid  flow  measurement  involves  several 
steps.  First,  it  is  necessary  to  "create" 
a  selected  plane  or  surface  within  the 
flow  field.  This  is  accomplished  by 
seeding  the  flow  with  small  tracer 
particles,  similarly  to  LDV  applications, 
and  illuminating  it  with  a  sheet  of  coher¬ 
ent  light,  as  shown  in  Figure  1.  A  pulsed 
laser  such  as  a  Ruby  or  a  NdYag  laser,  or 
a  CW  laser  with  a  shutter  is  normally  used 
as  the  light  source.  The  laser  sheet  is 
formed,  for  example,  by  focusing  the  laser 
beam  first  with  a  long  focal  length  spher¬ 
ical  lens,  to  obtain  minimum  thickness, 
and  then  diverging  the  beam  in  one  dimen¬ 
sion  with  a  cylindrical  lens.  The  light 
scattered  by  the  seeding  particles  in  the 
illuminated  plane  provides  a  moving 
pattern.  When  the  seeding  concentration 
is  low,  the  pattern  consists  of  resolved 
diffraction  limited  images  of  the  partic¬ 
les.  When  their  concentration  increases, 
the  images  overlap  and  interfere  to  pro¬ 
duce  a  random  speckle  pattern.  A  multiple 
exposure  photograph  records  this  moving 
pattern.  The  lower  particle  concentration 
originates  a  mode  of  operation  of  the 
technique  referred  to  as  Particle  Image 
Displacement  Velocimetry,  reserving  the 
term  Laser  Speckle  Velocimetry  for  the 
high  particle  concen t r a t i on  levels  where  a 
random  speckle  pattern  is  actually  formed 
(reference  6).  In  a  second  step  the  local 
fluid  velocity  is  derived  from  the  ratio 
of  the  measured  spacing  between  the  images 
of  the  same  tracer,  or  speckle  grain,  and 
the  time  between  exposures. 

Several  methods  exist  to  convert  the 
information  contained  in  the  multiple- 
exposed  photograph,  or  specklegram,  to 
flow  field  data  such  as  velocity  or  vor¬ 
ticity.  The  recorded  image,  whether 
formed  by  isolated  disks,  in  the  case  of 
low  particle  concentration,  or  speckle 
grains  for  high  particle  concen t ra t i on  is 
a  complicated  random  pattern.  It  would  be 
very  difficult  to  measure  the  local  dis¬ 
placements  by  visual  or  computer-aided 
inspection.  However,  it  is  important  to 
cealize  that  the  multiple  exposure  photo¬ 
graph  results  in  a  periodic  random  image 
from  which  the  periodicity  information  can 


he  retrieved  using  Fourier  or  Auto- 
coirelation  analysis.  basically,  the 
multiple-exposed  photographs  or  speckle- 
grams  can  be  analyzed  either  on  a  point- 
by-point  basis,  which  yields  measurements 
of  the  local  displacements  (velocity), 
(refs.  7-ll)  o  l  w  i  t  h  a  whole  field  filter¬ 
ing  technique,  which  yields  isovelocity 
contours  (ret.  9  )  .  Tile  method,  which  has 
been  selected  and  implemented  by  the  Fluid 
Mechanics  Research  Laboratory  at  the 
Florida  State  University,  is  the  Young's 
fringes  method.  The  local  displacement  is 
determined  using  an  focused  laser  beam  to 
interrogate  a  small  area  of  the  multiple 
exposed  photograph  transparency.  The  dif¬ 
fraction  produced  by  coherent  illumination 
of  the  multiple  images  in  the  negative 
generates  Young's  fringes,  in  the  Fourier 
plane  of  a  lens,  provided  that  the 
particle  images  correlate.  This  is  shown 
schematically  in  Figure  2.  These  fringes 
have  an  orientation  which  is  perpendicular 
to  tile  direction  of  the  local  displacement 
and  a  spacing  inversely  proportional  to 
the  displacement.  The  use  of  Young's 
fringes  eliminates  the  difficulties  of 
finding  the  individual  image  pairs  in  the 
photograph.  The  basis  of  the  Young's 
fringe  method  is  described  in  reference 
8  . 

The  photographic  recording  method 
discussed  above  lias  the  disadvantage  that 
the  photograph  consisLs  of  particle  pairs 
which  have  a  180  degrees  ambiguity  in  the 
direction  of  the  velocity  vector.  In 
addition,  it  has  been  shown  (reference  10) 
that  the  velocity  dynamic  range  of  the 
technique  is  limited  to  a  maximum  value  of 
about  10.  In  most  flows  of  interest  (e.g. 
boundary  layers  and  separated  flows),  this 
dynamic  range  is  not  sufficient  to  capture 
the  flow  field  in  its  entirety.  These 
limitations  are  critical  when  measuring 
complex  flows  having  flow  .  reversals  and 
s t agna tion  areas. 


A  method  to  resolve  both  the  ambiguity 
of  the  velocity  vector  as  well  as  to  im¬ 
prove  the  technique's  velocity  range  is 
incorporated  in  this  experiment.  This 
method  proposed  by  Lourenco1 1  and 
Adrian12,  commonly  known  as  "velocity  bias 
technique".  Consists  of  recording  the  flow 


experiment. 
Lourenco1 1 


range  is 
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field  in  a  moving  reference  frame,  thus 
superposing  a  known  velocity  bias  to  the 
actual  flow  velocity.  This  effect  may  be 
accomplished  in  several  ways,  in  particu¬ 
lar,  using  a  moving  camera  during  the 
photographic  recording  or  by  optical  means 
using  scanning  or  rotating  mirrors.  For 
the  data  presented  here,  a  rotating  mirror 
was  used  to  displace  the  image  during  the 
exposure  witti  a  pre-determined  velocity. 
A  schematic  of  the  rotating  mirror 
arrangement  is  shown  in  figure  3. 
Consider  two  particle  pairs  Ag Bg  and  Cg D0 
having  equal  displacements  in  opposite 
directions  in  the  object  plane.  By  intro¬ 
ducing  a  45°  mirror  between  the  camera 
lens  and  the  object  plane,  the  correspond¬ 
ing  displacements  appear  in  the  film  plane 
as  AB  and  CD  with  equal  magnitudes.  When 
the  mirror  is  rotated  by  an  angle  of  A  9 
between  exposures,  the  displacements 
corresponding  to  A  B  and  C  D  appear  in 
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the  film  plane  as  AD1  and  CD1  with 
different  magnitudes.  The  correct  dis¬ 
placement  or  velocity  with  its  direction 
can  now  be  obtained  upon  removal  of  the 
velocity  bias.  An  example  of  the  flow 
field  obtained  with  and  without  the 
velocity  bias  can  be  seen  in  figure  4. 
This  flow  represent  a  typical  separated 
flow  field  containing  flow  reversal  and 
stagnation  areas  (a  discussion  of  this 
flow  field  is  given  later). 

3.  Apparatus,  Instrumentation  and 
—  Procedures 

The  time-space  development  of  the 
unsteady  separated  flow  generated  by  an 
high  angle  of  attack  (oc=30°)  NACA  0012 
airfoil  impulsively  started  from  rest  is 
examined  using  Particle  Displacement 
Velocimetry.  The  flow  is  created  by 
towing  the  airfoil  in  the  reduced  scale 
Fluid  Mechanics  Research  Laboratory  towing 
tank  facility.  The  tank  is  300  x  200  x 
600  mm.  A  detailed  examination  showed 
that  the  motion  of  the  towing  carriage  is 
smooth  and  vibration  free.  The  towing 
carriage  is  driven  by  a  variable  D.C. 
motor,  and  the  towing  velocity  can  vary 
from  0  to  100  mm/sec.  For  the  photo¬ 
graphy,  a  35mm  camera  (Nikon  F-3)  is  used. 
To  photograph  the  flow  at  regular  time 
intervals,  the  photographic  camera  has  a 
electric  winding  device.  The  photographic 
time  interval  available  with  this  camera 
can  be  continuously  varied  up  to  a  maximum 
of  6  frames  per  second.  Two  options  are 
available  to  fix  the  camera;  one  by 
attaching  it  to  the  towing  carriage, 
which  means  an  observation  point  fixed  in 
relation  to  the  model,  and  the  other  by 
attaching  it  to  the  frame  of  the  water 
tank,  which  means  an  observation  point 
fixed  in  relation  to  the  fluid.  The 
selection  of  these  two  depends  upon  the 
flow  field  being  photographed. 

In  this  experiment  the  airfoil  chord 
is  60mm  and  is  towed  with  a  velocity  of 
22mm/sec.  The  corresponding  Reynolds 
number  was  1400.  The  fluid  used  in  this 
experiment  was  water  seeded  with  4,»m 
metallic  coated  particles  (TSI  model 
10087).  For  the  illumination,  a  laser 
beam  from  a  5  Watt  Argon-Ion  Laser 
( Spec t r a- Phy s i cs  series  2000)  is  steered 
and  focused  to  a  diameter  of  ,3mm  using  an 
inverse  telescope  lens  arrangement.  A 
cylindrical  lens,  with  a  focal  length  of 
-o.34mm,  is  used  to  diverge  the  focused 
beam  in  one  dimension,  creating  a  light 
sheet.  The  laser  sheet  is  70mm  wide  and 
illuminates  the  mid-span  section  of  the 
airfoil.  For  the  multiple  exposure,  the 
CW  laser  beam  is  modulated  using  a  Bragg 
cell.  In  this  experiment,  the  laser  power- 
density,  IQ ,  of  the  sheet  was  .27  W/mm2  . 
In  order  to  record  the  time  development  of 
the  flow  field,  the  camera  was  attached  to 
the  towing  carriage  and  the  frequency  of 
which  the  multiple  exposures  were  taken 
was  set  at  2.0Hz.  The  aperture  of  the 
lens  with  a  focal  length  of  50mm  and  a 
spacer  of  12mm,  was  set  at  F45.6  and  the 
resulting  magnification  factor  was  0.40. 
For  the  multiple  exposure,  the  time 
between  exposures,  T,  and  the  exposure 


time,  t,  are  chosen  according  to  the  cri¬ 
teria  discussed  in  reference  6.  The  time 
between  exposures  was  10msec.  For  optimum 
exposure,  file  exposure  time  was  lmsec, 
which  corresponded  roughly  to  (d,/MVmax), 
where  D  is  Lite  analyzing  beam  diameter,  V 
is  the  maximum  expected  velocity  plus  the 
shift  velocity  in  the  field  and  d  is  the 
particle  image  diameter  expressed  in  terms 
of 

=  (d2  +  d2)2 

with  d  the  particle  diameter  and  d  the 
edge  spread  caused  by  the  limited  response 
of  the  recording  optics  (ref.  6). 

Data  Process i ng 

The  fringe  images  were  acquired  and 
analyzed  using  the  digital  image  analysis 
system  of  the  Florida  State  University 
FMRL  (Fig.  5).  This  system  consists  of 
the  following  components:  a  DEC  LSI-11/73 
host  computer,  Gould  IP-8500  Digital  image 
processor  which  includes  four  memory  tiles 
for  storage  of  image  data  in  a  512  x  512 
format  with  a  resolution  of  8  bit  pet 
pixel,  a  frame  digitizer,  a  pipeline  pro¬ 
cessor  and  a  video  output  controller  to 
convert  digital  to  analog  information  for 
display  on  a  color  monitor.  The  system 
also  includes  a  two-dimensional  Klinger 
traversing  mechanism  with  a  controller  for 
the  purpose  of  automatically  scanning  the 
film  transparencies.  Two  methods  are 
available  and  used  for  fringe  analysis 
(ref.  10).  The  first  one  is  an  inter¬ 
active  method  in  the  sense  that  it 
requires  the  assistance  of  an  operator. 

The  inconvenience  of  the  one¬ 
dimensional  averaging  method  is  the  need 
for  an  external  adjustment  of  the  angle  of 
the  fringes  by  an  operator.  This  problem 
can  be  by  passed  by  computing  the  velocity 
components  along  independent  directions. 
Because  each  line  of  the  fringe  frame  can 
be  considered  as  a  noisy  periodic  signal 
with  variable  phase,  the  automatic  deter¬ 
mination  of  a  velocity  component  can  be 
performed  only  by  averaging  over  a  quanti¬ 
ty  independent  of  the  phase.  The  autocor¬ 
relation  for  each  line  or  its  Fourier 
transform  for  the  power  spectrum  satisfies 
ttiis  requirement.  The  m  velocity  compo¬ 
nent  can  be  computed  from; 

511  [E  ( I ( m, n ) I ( m+u , n )  ) 
g(u)  -  E  m _ _ -511<U<511 

n“°  E  (  I  (  m ,  n  )  ) 2 


This  algorithm  has  been  implemented 
using  the  pipeline  processor  of  the  Gould 
IP-8500  image  processor  to  perform  simul¬ 
taneously  the  autocorrelation  for  all  the 
lines  of  a  frame.  For  an  accurate  esti¬ 
mate  of  the  velocity  magnitude  and  direc¬ 
tions,  four  of  such  full  image  operations, 
yielding  four  autocorrelation  functions, 
are  required.  From  these  the  velocity 
vector  is  determined  by  selecting  the 
values  of  the  components  which  have  been 
comnuted  from  autocorrelations  having  the 
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highest  SNR,  and  visibility.  (The  compu¬ 
tation,  which  includes  the  de  te  irmi  na  t  i  on 
of  the  fringe  angle  and  position  updating 
of  the  film  transparency  scanning  mecha¬ 
nism,  is  completed  in  a  few  seconds, 
typically  4-5  sec,  using  the  PDP  11-71 
compute i . ) 

The  overall  accuracy  of  the  technique 
was  evaluated  using  a  method  described  in 
reference  10.  A  uniform  flow  field  is 
created  by  producing  a  multiple  exposure 
photograph  of  the  still  seeded  water,  in 
the  water  tank,  with  a  camera  moving  at 
constant  speed.  For  the  multiple  exposure 
photograph  a  number  of  time  between 
exposures  are  used,  thus  resulting  in 
photographs  with  particle  pairs  at 
different  known  distances. 

In  the  absence  of  a  systematic  bias, 
the  standard  deviation  of  the  obtained 
velocity  distribution  is  an  estimate  of 
the  mean  measurement  error.  From  the 
error  analysis,  it  is  believed  that  the 
velocity  data  is  obtained  with  an  accuracy 
of  2  percent  or  better. 

4,  Results  and  Discussion 

Typical  multiple  exposure  photographs 
of  the  flow  generated  by  the  impulsively 
started  NACA  0012  from  test  for  different 
times  are  shown  in  figure  6.  The  photo¬ 
graphic  arrangement  was  purposely  adjusted 
to  enhance  the  view  of  the  flow  field  on 
the  upper  surface  of  the  airfoil  rather 
than  to  show  the  entire  flow  around  the 
airfoil.  Consequently,  the  details  of  the 
flow  under  the  airfoil  can  not  be  seen 
clearly  in  these  photographs.  The  angle 
of  attack  of  the  airfoil  is  set  at  30°. 
These  pictures  display  the  flow  field  from 
the  leading  edge  to  a  downstream  location 
of  about  1.5  chords.  Photographs  were 
taken  at  a  frequency  of  2Hz  .  A  total  of 
34  pictures  were  obtained  covering  the 
range  of  t'  from  0.1  to  5.  The  non- 
dimensional  time,  t‘  (=  Ut/c,  where  U  is 
the  free  stream  velocity,  t  is  the  time 
between  two  successive  pictures,  and  c  is 
the  chord  of  the  airfoil)  between  suc¬ 
cessive  pictures  was  0.167.  However,  in 
figure  6,  only  a  limited  number  of  pic¬ 
tures  are  included.  The  quadruple  exposed 
photographs  shown  here  increase  the  SNR 
(signal  to  noise  ratio),  the  fringe  visi¬ 
bility,  and  provide  an  excellent  flow 
visualization.  The  added  advantage  of 
providing  a  good  flow  visualization  is  an 
asset  of  the  PIDV  technique. 

From  the  flow  visualization  pictures 
shown  in  figure  6,  the  following  obser¬ 
vations  are  made.  At  the  start  of  the 
aicfoil  a  vortex,  at  the  trailing  edge, 
commonly  known  as  "starting  vortex",  is 
generated  and  is  carried  away  from  the 
body.  concomitant  with  this  is  the 
generation  of  a  separation  bubble  at  the 
leading  edge  of  the  airfoil  (figure  4a). 
At  a  later  time,  for  example  at  t'  =  1.2 
(figure  4b),  the  separation  bubble  grows 
into  an  isolated  primary  vortex  with 
"secondary  vortices"  following  behind  it. 
Similar  type  of  vortex  structure  was  also 
observed10 in  the  flow  behind  a  circular 


cylinder.  This  multiple  vortex  structure 
continue  to  grow  together  until  the  t‘ 
reaches  a  value  of  about  2.5  (figure  4c 
and  4d).  At  t'  =  2.5  (figure  4d),  because 
of  the  close  proximity  of  the  primary 
vortex,  a  trailing  edge  vortex  is 
generated.  At  l'  =  2.75  the  primaiy 
vortex  abruptly  moves  away  from  the  sur¬ 
face  of  tile  airfoil  leaving  behind  a 
“vortex  sheet"  like  structure  (figure  4e). 
For  t‘  >  3.0,  this  "vortex  sheet"  rolls  up 
into  distinct  vortices  and  they  grow  in 
size  with  time  as  shown  in  figure  4f  -  4g. 
During  this  process,  the  trailing  edge 
vortex  also  grows  and  as  a  result  the 
whole  flow  field  becomes  very  complex. 
Close  to  the  surface  of  the  airfoil,  a 
small  vortex  remain  present  for  t‘  >  3.0. 
This  vortex  has  the  same  sign  of  rotation 
as  the  trailing  edge  vortex.  A  similar 
vortex  structure  was  observed  by  Ho1 J ,  who 
calls  it  an  "induced  vortex"  and  associ¬ 
ates  it  with  unsteady  separation 
phenomenon . 

The  velocity  data  is  acquired  in  a 
Cartesian  mesh  by  digital  processing  of 
the  Young's  fringes,  produced  by  point-by¬ 
point  scanning  of  the  positive  contact 
copy  of  the  photograph.  The  scanning  step 
size  and  the  dimension  of  the  analyzing 
beam  are  0.5mm,  which  corresponds  to  a 
spatial  resolution  of  about  1.25mm  in  the 
object  plane  or  about  0.02c.  The  fringes 
were  processed  using  the  method  described 
in  the  previous.  section.  The  resultant 
two-dimensional  velocity  fields,  corre¬ 
sponding  to  figure  6a  -  6h,  are  shown  in 
figure  7.  The  length  of  each  vector  is 
proportional  to  the  local  velocity  at  that 
point.  Tire  color  code  superimposed  on  the 
velocity  data  represents  the  vorticity 
level,  the  magnitude  of  which  is  given  by 
the  color  bar  at  the  bottom  of  the 
picture.  The  red  and  green  colors  repre¬ 
sent  the  peak  positive  and  negative  vor¬ 
ticity  regions  respectively.  This  type  of 
display  clearly  depict  the  various  regions 
of  vorticity  and  its  strength.  As  dis¬ 
cussed  above,  the  presence  of  primary 
vortex,  secondary  vortices,  vortex  sheet 
trailing  edge  vortex  and  the  induced 
vortex  on  the  surface  of  the  airfoil  are 
clearly  depicted  along  with  their  time- 
space  development.  The  detailed  analysis 
of  this  data  is  being  conducted  at  this 
time  and  will  be  reported  later. 

5.  Conclusions 

A  recently  developed  velocity  measure¬ 
ment  technique,  known  as  Particle  Image 
Displacement  Velocimetry  (PIDV),  has  been 
briefly  described.  Using  this  technique, 
the  time-space  evolution  of  the  flow 
generated  by  an  impulsively  started  high 
angle  of  attack  (a=30°)  NACA  0012  airfoil 
is  presented.  This  experiment  illustrated 
the  technique’s  capabilities  to  record 
with  accuracy  the  complex  unsteady 
separated  flow. 

The  technique  has  been  shown  to  pro¬ 
vide  both  flow  visualization  and  quanti¬ 
tative  measurements,  which  include  the 
velocity  and  vorticity  fields. 
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The  unsteady  separated  flow  field 

generated  by  an  high  angle  of  attack  air¬ 
foil  contains  many  large  scale  vortical 
structures  such  as;  primary  vortex 

generated  at  the  leading  edge  with  second¬ 
ary  vortices  upstream  of  it,  trailing  edge 
vortex,  vortex  sheet  and  an  induced  vortex 
in  the  upper  surface  of  the  airfoil.  The 
origins  and  time  development  of  these  are 
clearly  depicted  by  the  instantaneous 

velocity  and  vorticity  fields  obtained 
using  PIDV. 
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Figure  3.  Rotating  Mirror  arrangment  for  the  Velocity  Bias  Technique 
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Figure  4.  Instantaneous  Velocity  Field  a)  Before  removal  of  velocity 
bias  b)  after  removal  of  the  velocity  bias. 


